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formed this research for the Earth Resources Research Division, Space and Life 
Sciences directorate. National Aeronautics and Space Administration, at the 
Lyndon B. Johnson Space Center, Houston, Texas. Mrs. Kathy McIntyre provided 
the programming support for accessing the merged data and ground- truth label 
files on tape. 
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1. INTRODUCTION 


Recently, considerable Interest has been shown In developing techniques 
(refs. 1 and 2) for the analysis of multichannel Imagery data (such as remotely 
sensed mul t1 spectral scanner data acquired by the Landsat series of satellite) 
for Inventorying natural resources, predicting crop yields, detecting mineral 
and oil deposits, etc. One of the Important objectives In the analysis of 
remotely sensed Imagery data Is to estimate the proportion of the crop of 
Interest In the Image. Nonsupervi sed classification or clustering techniques 
(ref. 2) which partition the Image Into Its Inherent modes or clusters have 
been found to be effective In the classification of Imagery data for proportion 
estimation. 

Usually, agricultural Imagery data have a fleld-llke structure (ref. 3). The 
resolution element or pixel of the remote sensing Imagery corresponds to 
approximately 0.44 hectares (l.I acres) on the ground. A significant portion 
of the Imagery data will contain mixture pixels (I.e., pixels containing 
objects from more than one class) whenever the objects being viewed by multi* 
spectral scanner (MSS) are not large enough relative to the size of a resolu- 
tion element. The percentage of mixture pixels In the Image depends In general 
on the size of the fields. By analyzing a nunber of remotely sensed multi- 
spectral agricultural Images, Nalepica and Hyde (ref. 4) have estimated that, 
for 20-acre fields, the percentage of mixture pixels In the Image Is around 
40 percent i and, for fields between 60 acres and 100 acres, the percentage of 
mixture pixels exceeded 20 percent. Hence, to be able to accurately estimate 
the proportion of the crop of Interest In the Image, It Is necessary to deal 
with the mixture pixels. 

Recently, several researchers (refs. 5 and 6) have attempted to partition or 
segment a multichannel Image Into pure pixel (I.e., pixels containing objects 
of a single class) regions or fields and Into mixed pixel or boundary pixel 
regions. There Is considerable Interest In developing techniques (refs. 4 and 
7) for estimating the proportion of classes In the mixed pixels. In all the 
proposed methods the proportions of classes In the mixed pixels are estimated 
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as follows. Assuming the spectral response vector of the mixed pixel as 
Gaussian, the proportions of classes In the mixed pixel are estimated as those 
that maximize the likelihood of occurrence of Its spectral response vector. 

One of the reasons these approaches are not successful. In general. Is that the 
Individual observation vectors are noisy. In this paper, techniques are 
developed for estimating the proportions In the mixed pixels by the characteri- 
zation of region of mixed pixels. The probability density function of the pro- 
portion of classes In the mixed pixels Is estimated using Information from the 
spectral vectors of a set of mixed pixels from the mixed pixel region. 

Estimates for the proportion of classes In the mixed pixels are then obtained. 

This paper Is organized In the following manner. In section 2, relationships 
are developed between the moments of the spectral vectors of the mixed pixels 
and the moments of the spectral vectors If the entire pixels contain objects 
from a single class. In section 3, expressions are developed for the maximum 
likelihood estimates of the parameters of the probability density function of 
proportions of classes In the mixed pixels of the region when the mixed pixels 
contain two classes of objects. Experimental results from the processing of 
remotely sensed agricultural Imagery data are presented In section 4. In 
section 5, the results of section 3 are generalized when the mixed pixels 
contain more than two classes of objects. The concluding simmary Is given In 
section 6. In the appendix A, using Information from the spectral vectors and 
the proportion of classes of a given set of mixed pixels from the region, 
expressions for the estimation of the probability density function of the 
proportion of classes In the mixed pixels of a region are presented. The 
dependencies between the spectral response vectors of the subpixels of the 
classes are dealt with In appendix B. 
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2. A MODEL FOR THE CHARACTERIZATION OF BOUNDARY PIXEL REGIONS 


A model for the characterization of regions of boundary pixels Is presented In 
this section. With the present resolution, a pixel or a resolution element of 
Landsat MSS Images represents an acre of ground area. It Is observed that more 
than 30 percent of the pixels of a typical MSS Image are boundary or mixed 
pixels (I.e., pixels containing more than one class of objects). 

Let a pixel consist of K small cells of equal size, and let be the number of 
cells containing the 1th class. Let x^j be a random vector representing the 
spectral response of class 1 In the jth subcell of these cells. The situa- 
tion Is Illustrated In the figure 2-1, where for convenience the subcells of 
class 1 are shown as a contiguous block. 


Class 1 cells 


Ml 

M2 





Mj 





Mk^ 


Figure 2-1.- Spectral response vectors associated with 
the cells of class 1 In a resolution element. 


Let the spectral response vectors j ■ 1, 2, ••*, K^, have mean M{ and 
covariance matrix for 1 • I, 2, •••, R, where R Is the nunber of classes of 
objects In the resolution element. Let the total response for the resolution 
element be represented by the random vector X. Assume that X can be written as 


X ■ 



'U 


( 2 - 1 ) 
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Let K be the total nueiber of subcells of the resolution element, where 


K « 



( 2 - 2 ) 


If the entire resolution element were to consist of class 1, assuming 
Independence between the spectral response vectors of the subcells,^ the mean 
vector and the covariance matrix of X can be obtained as follows* 


and 


« E(X) - KM) 

• COv(X) » KE) 


(2-3) 


Since there are actually K^ subcells of the 1th class, the mean of X Is 

£(X) « ^ K^M) ■ ^ a^KM) 

» ^ » M(o) (2-4) 

1-1 ' ^ 


where 


*1 “IT 


(2-5) 


and Is the proportion of class 1 In the resolution element. The proportions 
satisfy the following relationships. 


and 


a. > 0 



1 - 1 . 2 . 



( 2 - 6 ) 


^The dependencies between the spectral response vectors of the subpixels of the 
classes are dealt with In appendix B. 
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If th« random vectors associated with the subcells of different classes are 
also assumed to be Independent* the covarUnce matrix of X can be written as 


cov(X) « 



E(a) 


* (2-7) 


Let the elements 1 « 1* 2* •**, R* of the vector a satisfy equation (2*6). 
Let p(o) be the probability density function of a characterizinq a reqlon of 
mixed pixels* Let be the reqlon of a In which the constraints of equation 
(2-6) are satisfied. Let p^(X) be the probability density function of the 
spectral response vectors X of the mixed pixels. It can be written as 

P,(X) . ^ P,(X..)<1. 

01 

• f p,_(X| a)p( a)da (2-8) 

.a 

a 


One of the Important objectives In the analysis of remotely sensed Imagery data 
Is to estimate the proportion of the class of Interest In the Image. If p(a) 

Is known or estimated, given an observation vector X of a mixed pixel, the 
Bayes posteriori estimate for the proportion of classes In the mixed pixel can 
be obtained as follows. 


a > E(alX) 


<q)„( a|X)via 
In 

a 


f flp_(X|o)p(o)do 

•'a " 

tt 


/ P_(X|a)p(a)da 

rt In 



(2-9) 


2-3 


3. ESTIMATION OF p(a) VIHEN THE MIXED PIXELS CONTAIN TWO CLASSES OF 09JECTS 


The problem of estimation of p(a) to characterize a region of mixed pixels, 
given the spectral response vectors of a set of mixed olxels from the region. 

Is considered In this section. Very often the proportion of classes In the 
mixed pixels Is unknown. The Identification of mixed or border pixels, 
however, can be obtained by using either the clustering algorithms or the 
segmentation algorithms. Assuming functional forms for p(a), expressions are 
developed In the following paragraphs for obtaining the maximum likelihood 
estimates of the parameters of p(a) using Information from the observation 
vectors of a set of mixed pixels. From the analysis of several ground-truth 
Images, It Is observed that suitable functional forms for p(a) are (a) the beta 
distribution function and (b> the density function representing the portion of 
a Gaussian curve In the region of Interest. These functional forms are 
described In the following paragraphs. Let a be the proportion of class 1 In 
the mixed pixel. Then (1 - a) Is the proportion of class 2* 


a. Beta distribution: Modeling p(a) as a beta distribution In terns of 
unknown parameters. It can be written as 


p(ol * 


Aa’^d-al^^ ; 0 < o < 1 
0 ; elsewhere 


where b > -1 and c > -1 are the parameters to be estimated and the 
constant A Is given by 

. _ r(b ♦ c ♦ 2) 

A • TTFTTTrrrnT 


(3-1) 


(3-2) 


and r(*) is a usual gamma function. 


b. Gaussian surface: The probability density function p(a) can also be 
modeled as a portion of Gaussian surface In the illowable region of a. 
That Is, p(a) can be written as 


p(al 


JT^ fU) dc 
0 ; otherwl se 


(3-3) 
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whert f(a) It a Gauttlan density function with man ay and variance S^. The 
paramters aif and are to be estimated* The probability density function 
p(a) Is Illustrated In figure 3-1* 



Figure 3*1. - The probability density function p(a), when 
modeled as the portion of a Gaussian surface. 

3.1 MAXIMUM IIKEUHOOO ESTIMATION OF p(a) 

The estimation of the paramters of p(a) Is formulated In this section as that 
of a mximum likelihood estimtion problem. It Is assumd that the spectral 
response vectors t • 1, 2, N, of a set of mixed pixels are given. The 
log likelihood of the occurrence of the set of given obsenatlon vectors can be 
written as follows. 


L • ^ logCpjXj] 

• Pj,(X^|a)p(a)daj (3-4) 

Closed form solutions for the paramters of p(a) that mximize L seen to be 
difficult when the functional form of equation (3>1) or equation (3-3) Is used 
for p(a). In general, the paramters of p(o) that mximize L can be obtained 
using optimization techniques such as the Oavldon-Tletcher-Powell procedure 
(refs. 8 and 9). However, Iterative equations, which are similar to mximum 
likelihood equations In clustering (refs. 10 and 11), for the estimtion of 
paramters of p(a), can be obtained using the functional form for p(a) given by 
equation (3-3). The following maximum likelihood equations can easily be 
derived by differentiating L with respect to the parameters of p(o) and 
equating the resulting expressions to zerc. 
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•"f 



«q)|^(X^ |«)f(o)d« 

0 


(m- - a)f(a)da 


0 ^ 

p^(X^|a)f(a)do 

.^0 


/‘f(«)da 
. ^0 


(3-5) 


and 


Sf 




(a - P|^(X^ |a)f(a)da 

pjX^la)f(a)da 


♦ 


j[ [^f * 

fU)da 


f(a)da 


(3-6) 


Tbc use of equations (3-5) and (3-6) requires* In general* that the Integration 
be performed numerically. From equations (2-4) and (2-7)* for a particular a* 
the mean and the covariance matrix of the spectral vectors of the mixed pixels 
are given by the following. 

M(a) ■ oMj ♦ (I - a)M 2 (3-7) 
t(a) • o£^ ♦ (1 - a):^ (3-8) 


The resolution elements that contain a single class are called pure pixels. In 
the following equations* It Is assumed that the spectral vectors of the pure 
pixels are Gaussian. For a given a* It Is also assumed that the spectral 
vectors of the mixed pixels are Gaussian. In the estimation of ay end S^* by 
Iteratively using equations (3-5) and (3-6)* the computation can be consider- 
ably reduced by transforming the spectral vectors with a transformation matrix 
that simultaneously diagonalizes the covariance matrices Cj and £ 2 * Let A be 
the transformation matrix. Then we have [123 


Mtiere fJ • The matrices a and a are the eigenvalue and eigenvector 

matrices of The matrices A and a are the eigenvalue and eigenvector 
matrices of K, where 

K « (3-10) 

Let the spectral vectors be transformed Into vectors Y^, where 

« AX^ ; A • 1, 2. M (3-11) 

Let the means of the pattern classes be transformed Into where 

« AM^ ; 1 « 1, 2 (3-12) 

From equations (3-7), (3-8), (3-9), and (3-12), for a given o, the mean and the 
covariance matrix of the transformed spectra) vectors of the mixed pixels are 
given by the following. 

u(a) “ auj + (1 - a)ii2 (3-13) 

and S(o) * oI + (1 - o)A (3-14) 

The use of Pn|(Y{|a) In equations (3-5) and (3-6) reduces the computation con- 
siderably since the determinant and the Inverse of matrix S(a) can be computed 
directly from equation (3-14). An estimate for the proportion of the class of 
Interest (say class 1) In a mixed pixel with the transformed observation vector 
Y Is given by the following. 


a 


ap^(Y|a)f(a)do 
0 

p^(Y|a)f(a)da 


(3-15) 


3.2 MAXIMUM LIKELIHOOD ESTIMATION OF p(a), WITH THE CRITERION OF A LOWER BOUND 

mr 

It Is observed that In equations (3-5) and (3-6) the ntinerical Integration Is 
to be performed at each Iteration for the transformed spectral vector of every 
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given mixed pixel. In the following paragraphs, it is shown that the computa- 
tion can be considerably simplified by using a lower bound on the likelihood 
function as a criterion. By noting that the logarithm is a convex upward 
function, a lower bound on L of equation (3-4) can be obtained as follows. 

L > 4 (3-16) 


where 



A( a)p( a) da 


(3-17) 


and 


A(a) 



logCp„(Y^lo)] 


(3-18) 


3.2.1 MAXIMUM LIKELIHOOD EQUATIONS FOR THE ESTIMATION OF PARAMETERS OF p(o) 


The maximum likelihood equations for the estimation of parameters of p(o) that 
maximize L^ of equation (3-17) can easily be shown to be the following. 


"f ‘71 * rl 

A(a)f(a)do Jq f(a)d 


(3-19) 


and 


Sf * 


(a - m^)^A(a)f(a)da [s^ - (a - m^)^]f(a)da 

•'ft 


7 


A(a)f(a)da 


T 


( 3 - 20 ) 


f( a)da 


It is seen that the use of equations (3-19) and (3-20) requires the integration 
to be performed numerically, once for every iteration. In the transformed 
space, an expression for A(a) is given by the following. 

A(a) - log(2») - J 1<>9 IS(a)| - .^rCS’^IolSV] 

♦ SMV^(a)u(o) - ^"^(alS'^ialuIa) 

• ^ J i log[o ♦ (1 - a)X^] + ^ + (1 - o)x/j 
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where 


'■(iN 

a, 

• <»*n - U2i)CSM(i) - Nu2^] 

Ci » U2^[sM(i) - 5t«2^] - 7 SVd.i) 


(3-22) 


The diagonal elements of the eigenvalue matrix A are and the dimensionality 
of the patterns Is n. 


3.2.2 CLOSED FORM EXPRESSIONS FOR THE INTEGRALS IN EQUATIONS (3-19) AND (3-20), 
WHEN THE COVARIANCE MATRICES OF THE CUSSES ARE EQUAL 


In the following paragraphs, expressions are derived for the computation of the 
Integrals In equations (3-19) and (3-20) when the covariance matrices of the 
classes are equal. If the covariance matrices of the classes are equal, then 
« 1 for all 1 and A(a) In equation (3-21) becomes 

2 

A(o) » aa + bo + c (3-23) 


where 


and 


c 


a « 


b » 


1«1 



(3-24) 
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Let 


**p(' 7^)^^ 

The following can now easily be derived* 

jT f(a)da ■ 



(3-25) 


(3-26) 




(3-27) 


X I , Sr I -W. /-">f 1 (1 - "If) 

(o - m^)^f(a)do * ^ exp 1 — — exp 


/2T I /57 \‘^f/ ^ 


(1 - m^)' 
' ■ 2!Tr ■ 


J f (a - fflr 

o ^ 


)^f(a)do 


s3/2(/ m2\ U\ f (l-n,r)2 

— ^12 ♦ J.lexol-,1 - 2+ — .r-l— 


exp 


+ jT f( o)da 
0 


(3-28) 
(1 - w^)^ 

(3-29) 




3(1 - nif) (1 - m^)' 




*~7J7T 


exp 


(1 - m^)‘ 
5K>r' 


I + sj f(o)do 


(3-30) 


The Integrals In equations (3-19) and (3-20) Involving the term A(a) can be 
expressed In terms of the above equations as follows. 
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A(a)f(a)do ■ a r (a - m.)^f(a)do * Z(mf + b) (a - m.)f(o)da 
® •'o ^ 

♦ |ain^ + bm^ ■*' ^) f(«)da 

oA(a)f(a)da • i ^ - ■f)^f(a)da ♦ (3 m^ ♦ b) (a - a^)^f(a)4a 

* * 2tMf * ' «y)f(a)da * 4- ca^j f(a)da 


( 3 - 31 ) 


( 3 - 32 ) 


J (o - m^)^A(a)f(o)do ■ a (a - m^)^f(a)da+ (2am^ + b) J" (a - m^)^f(o)da 

♦ |am^ ♦ bm^ + c| J* (a - m^)^f(a)da 


( 3 - 33 ) 
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4. EXPERIMENTAL RESULTS 


In this section, some results from the processing of remotely sensed MSS Imagery 
data are presented. Several segments were processed In the following manner. 
For every segment, several acquisitions are acquired and the Images are regis* 
tered. Each acquisition Is a 4-channe1 Image. The 4-channe1 Image values are 
transformed Into greenness and brightness space (ref. 13), thus generating a 
2>channe1 Image. Two classes are considered. Class 1 Is wheat and class 2 Is 
pasture. The class of Interest In the Image Is wheat. The resolution element 
or pixel of the Image corresponds to approximately an acre on the ground. Each 
pixel Is divided Into six subpixels, and the true class labels, or the ground- 
truth labels for each of the subpixels, are acquired. The pixels containing 
only wheat, the pixels containing only pasture, and the mixed pixels having 
wheat and pasture In different proportions are located In the segment. The 
spectral response vectors of pure pixels are assumed to be Gaussian. For a 
given a, the spectral response vectors of the mixed pixels are also assumed to 
be Gaussian. Assuming the functional form of equation (3-3) for p(a), the 
maximum likelihood estimators for the parameters of p(a) are obtained using 
equations (3-S) and (3-6). The spectral vectors are transformed using a trans- 
formation matrix that simultaneously diagonalizes the covariance matrices of 
the two classes. Simpson's rule Is used for computing the Integrals numeri- 
cally. The proportion of classes of Interest (I.e., wheat) In the mixed pixels 
Is estimated using equation (3-15). The number of pixels from each of the 
classes and the number of mixed pixels are listed In table 1. Also Included In 
table 1 Is the average true proportion of wheat In the mixed pixels estimated 
from the ground- truth labels of the subpixels of the mixed pixels. The esti- 
mated proportion of wheat In the mixed pixels using equations (3-5), (3-6), and 
(3-15), after first Iteration and after the convergence, are listed In table 1 
for n « 2 and 4. For a subset of the segments of table 1, the estimated pro- 
portion of wheat In the mixed pixels Is listed In table 2 for n « 6 and In 
table 3 for n » 8. In general. It Is observed that the better 



A segment Is a 9- by ll-k11ometer (5- by 6-naut1 cal -mile) area for which the 
MSS Image Is divided Into a 117-row by 196-column rectangular array of pixels. 
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TABLE 1.- ESTIMATES OF PROPORTION OF WHEAT IN MIXED PIXELS FOR n = 2 AND n 


5 

1 

i 

1 

5 

<iodoe5ddooo 





1 i ! 2 ! 5 i 1 ! s 

<M 

• 

• 

o 

M 

1 

» 

S 

• 

o 

s 

s 

■ 

C 

1 1 1 i 1 1 M M 

e'dddddoeee 

1 

» 

s 

d 

i 

m 

• 

o 


1 

i i s s i I 1 S 1 1 

^ d d d d d o e e e 

• 

s 

wm 

• 

O 

<si 

1 

IS 

1 

• 

o 

1 

i I M M M M 

OOOOOOOOOO 

*4 

! 

• 

O 

7 

• 

e 

! 

2 

& 

0 

1 

e 

3 

H 

s 

l2Sg228Saa 


31 

T 

3 

« 

a. 

§§§§!§§§§§ 

• 

1 


s 

*4 

3 

a 

a 

¥1 

1 

U 

il 2S t S 1- 5 M »i J M li 

Is II II II II Is Is II ll 


[ 

sssssssiii 

«*#*•*•• 

M 

Z 



4-2 














TABLE 2.- ESTIMATES OF PROPORTION OF WHEAT IN NIXED PIXELS FOR 
















TABLE 3.- ESTIMATES OF PROPORTION OF WHEAT IN MIXED PIXELS FOR 





















proportion estlnatos are obtained for n * 4. It Is thought that the degrade* 
tion In the estimates with the Increase In the number of acquisitions Is due to 
the registration errors. 
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5. ESTIMATION OF p(a) WHEN THE MIXED PIXELS 
CONTAIN MORE THAN THO CLASSES OF OBJECTS 


Th« problem of estimation of p(a) when the mixed pixels contain more than two 
classes of objects Is considered In this section. The functional forms that 
can be used for p(a) are the multivariate generalization of the ones presented 
In section 3. These are described In the following paragraphs. 

a. 'Hie Oirichlet Distribution: If p(a) can be represented as a Oirichlet 

distribution function. It can be written as 

R a. 
p(a) • K n «| 

1-1 ^ 


■ K 


(■ 




( 5 - 1 ) 


where 



1, > 0 for 1 « 1, 2, •••, R, and 


K 



n Cr(a. ♦ D] 

J-l ^ 

The set of parameters (a^) are such that a^ > -1 for 1 
are to be estimated. 


(5-2) 


1, 2, •••, R, and 


b. 


The multivariate Gaussian surface: By modeling p(a) with the surface of a 
multivariate normal distribution In the region p(a) can be written as 


p(a) 


f(tt) . 

iJ^ f(a)da 


If a e P 


0 ; otherwise 


(5-3) 


where f(a) Is a Gaussian density function with the mean vector M^ and the 
covariance matrix tf. The parameters M^ and are to be estimated. 
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5.1 HAXIMUM tIKELlHOOO gSTlMATION OF p(g) 


61v«n th« spectral response vectors i ■ 1, 2, •••, N, of a set of mixed 
pixels, the I09 likelihood of the occurrence of the given set of observation 
vectors can be written as follows. 


I • log 


[A 

■ ^ p^(Xi I a)p( a)daj 


(5-4) 


In general, using the functional forms for p(«) that are given either In equa- 
tion (5-1) or In equation (5-3), the parameters of p(o) that maximize L can be 
obtained using optimization techniques such as Oavldon-Fletcher-Powell (refs. 8 
and 9). If the functional form given by equation (5-3) Is used for p(a), the 
following maximum likelihood equations for the estimation of parameters of p(a) 
that maximize L can easily be derived. 


y* ap|i|(X^ |a)f(o)da 
^0 

4 

f (M- - c>)f(a) 4 o 
0 

J" Pj,(X^ la)f(a)da 


f f(a)da 
*'a 



« J 


and 




"'•'ll 


^ (a • N,l(« • N^)\(Xf|a)f(«)e«* 


(5-5) 


f |iy|(X||a)f(a)da 


• (a - Nf)^jf(a)a« 


/ f(a)4« 


(5-6) 


It Is noted that In equations (5-5) and (5-6) the Integrals need to be computed 
for every spectral vector at each Iteration. 


5-2 


5.2 WAXIMUM IIKEIIHOOO ESTIWATION OF p(a) WITH THE CWITERION OF A LOWER BOUMO 

Mr 


since the loyerlthn Is a convex upward function, a lower bound on L of 
equation (5-4) can be c^ta1ned as 


where 


and 


L > L, 


Lj ■ A(a)p(a)da 
*.(o) • ^ logCp„(x^|a)] 


(5-7) 

(5-8) 

(5-9) 


If the functional form given equation (5-3) Is used for p(a), the following 
maximum likelihood equations for the estimation of parameters of p(a) that 
maximize can easily be derived. 


f oA(a)f(a)da f (N« - a)f(o)da 

‘'a ^ 

^ ^ o 

^ A(a)f(a)do f( a)do 

a Cl 


(5-10) 


and 


(a - M^)(a - M^)"^A(o)f(a)da jj [l^ • ( o - H^) (a - M^)^]f(o)do 

Jq A(a)f(a)da f( a)da 

a a 


(5-11) 


It IS Observed that the use of equations (5-10) and (5-11) requires the 
Integrals to be computed once for every iteration. 
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6. CONCLUDING SUMMARY 


Ont of tNo Important objectives In the processInR of remotely sensed Imagery 
data Is the estimation of the proportion of the class of Interest In the Image. 
Agricultural Imagery data usually have a fleld-llke structure. From the analy- 
sis of several Images. It Is observed that depending on the field size. a. sig- 
nificant portion of the Image contains mixed pixels d.e.. pixels containing 
more than one class of objects). Techniques are currently being developed for 
partitioning an Image Into regions of pure pixels and mixed pixels. 

This paper addresses the problem of estimating the proportions of classes In 
the mixed pixels. Relationships are developed between the moments of spectral 
vectors of mixture pixels and the moments of spectral vectors of pure pixels of 
different classes of objects as a function of the proportion of clssses In the 
mixed pixel. The probability density function p(a) of the proportion of classes 
In the mixed pixels can be modeled either as a Oirichlet distribution or as a 
normalized Gaussian surface In the region of Interest In terms of unknown param- 
eters. Given the spectral vector of a mixed pixel, the proportion of classes In 
the pixel can then be estimated using p(a). 

By modeling p(a) as a normalized Gaussian surface, expressions are developed 
for obtaining the maximum likelihood estimates of Its parameters using Informa- 
tion from the spectral vectors of a set of mixed pixels from the mixed pixel 
region. These Involve evaluation of Integrals numerically. If the mixed pixels 
contain two classes of objects, the computation can be considerably reduced by 
transforming the spectral vectors using a transformation matrix that simultane- 
ously diagonalizes the covariance matrices of the classes. Closed form expres- 
sions are developed for the computation of the Integrals when the covariance 
matrices of the classes are equal. 

Experimental results are presented from the processing of remotely sensed 
agricultural Imagery data. Two classes are considered. Class 1 Is wheat and 
class 2 Is pasture. The mixture pixels contain wheat and pasture In different 
proportions. The proportion of wheat In the mixture pixels Is estimated using 
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equations (3-5), (3-6), and (3-15) with the dimensionality of the spectral 
vector n where n « 2, 4, 6, and 8. The proportion of wheat In the mixed pixels 
Is also estimated using the true class labels of the six subpixels of a pixel 
from the ground truth. Better proportion estimates are obtained for the value 
of n « 4 or with the data from two acquisitions. 

Expressions are developed for the estimation of parameters of p(a) based on the 
criterion of the minimum sum of the squares of errors, using Information from 
the spectral vectors and proportions of classes of a set of mixed pixels from 
the mixed pixel region. Furthermore, the effect of the dependencies between 
the spectral vectors of the subpixels of the classes on the moments of the 
spectral vectors of the mixed pixels 1$ considered In appendix B. 
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APPENDIX A 

ESTIMATION OF p(o) WITH THE CRITERION OF THE 
MINIMUM SUM OF THE SQUARES OF ERRORS 



APPENDIX A 

ESTIMATION OF p(o) WITH THE CRITERION OF THE 
MINIMUM SUM OF THE SQUARES OF ERRORS 

The problem of characterization of a region of mixed pixels through the est1ma> 
tion of p(a) using Information from the spectral vectors of a given set of 
pixels from the region was treated In sections 3 and 5. If the proportion of 
classes In the mixed pixels of the given set are also known, the problem of 
estimation of p(a) using all the available Information (In addition to the 
spectral vectors) Is considered In this appendix. The estimates of the param- 
eters of p(o) are obtained using the mlnlmtin sun of the squares of errors as a 
criterion. 

Let Xjj, t » 1, 2, N, be the n-dimenslonal spectral response vectors of the 
given set of rdxed pixels. Let a^, 1 » I, 2, N, be the R-dImenslonal 
vectors of proportions of classes In the mixed pixels. Given the spectral 
vector of a mixed pixel and p(a), an estimate for the proportions of classes 
In the pixel Is given by 

f ap(X Ja)p(a)da 

\ 

a. » —2 (A-1) 

f p(X. |a)p(a)da 


If the functional form of equation (5-3) Is used for the probability density 
function p(a), equation (A-1) can be written as 


f op(X. |a)f(a)da 

* 

f p(X|Ia)f(.)da 
a 


(A-2) 


where f(a) Is a Gaussian density function with the mean vector Mf and the 
covariance matrix £f. The criterion of the minimum sun of the squares of 
errors can be used for obtaining the parameters Mf and tf of p(a). The sum of 
the squares of errors, e, of the proportion estimates can be written as 
f ol 1 ows . 


A-1 


(A-3) 


' ■ |j (“i - M 

If 0 is 4 parameter of p(a), differentiating e with respect to 0 results In 



From equation (A-2). we get 


3a^ 

W 


J ap(X^ |a) [f(a)]da - J p(X^ I a) ^ [f( o)]dc 

o 

p(X^|a)f(g)da 


(A-5) 


(A-6) 


Differentiating f(a) with respect to Its mean vector Mf yields 

tj^g - M^)f(g) 

Let be the elements of the matrix Differentiating f(g) with respect 

to results In the following. 


3f( 

3v 


— • “it - *“t T- 


(A-7) 


and 




Where o^j are the elements of the matrix if and Is the 1th element of the 
vector Mf. Substitution of equations (A-5), (A-6), and (A-7) In equation (A-4) 
yields Iterative equations (A-8) and (A-9), which are similar to maximum likeli- 
hood equations (refs. ID and U), for the estimation of parameters Mf and if of 

p(a). 


When there are only two classes In the mixed pixel, as shown In section 3, the 
computation can be greatly reduced by transforming the spectral vectors using a 
transformation matrix that simultaneously diagonalizes the covariance matrices 
of the classes. 
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a p(X . I a)f ( a)da 


00 
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APPENDIX B 

EFFECT OF CORRELATIONS BETWEEN THE SPECTRAL VECTORS OF SUBPIXELS 
ON THE MOMENTS OF SPECTRAL VECTORS OF MIXED PIXELS 


In section 2, It Is assumed that the spectral vectors of the subpixels are 
Independent. The purpose of this appendix Is to take Into account the correla- 
tions between the spectral vectors of subpixels In developing expressions for 
the moments of the spectral vectors of the mixture pixels. If the entire 
resolution element were to consist of class 1, the spectral vector X of the 
resolution element can be written In terms of the spectral vectors of the 
subpixels as 


X 



*1J 


(B-1) 


The mean vector M^ and the covariance matrix of X can be obtained as 
follows. 


M^ » E(X) » KMj 


(B-2) 


» cov(X) 


E^ (X.. - Mj)(x.. - Ml)^ + Z (x.. - Ml)(x.^ - m;)^ 
J-1 ^ ^ j«l k*l ^ IK 1 

L 

kt»J 


(B-3) 


If the spectral vectors of the subpixels are Independent, the second term on 
the right-hand side of equation (B-3) becomes zero. Let 



Let the covariance matrix of the random vector Mhich can be 

written as 


Let 




‘•Isr 


E' 

^Isr 




^Isr 


Msr 
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(B-6) 


If the random vectors x^^ and x^ ^ are Gaussian with mean M| and covariance 
matrix £|, the conditional probability density p(x^^|x^^) Is normal with mean 

vector M| - - M^) and covariance matrix Now consider 

'['■'ir - - "i>^] ■ /<”lr - "!> [/'"ts - 

(B-7) 


Using equation (B-7) In equation (B-3) yields 


Kc: 


-'i 1 1 

r»l s»l 


s-1 

srr 


^Isr^ls 


(B-8) 


It Is assumed that the covariance matrix of radiance vectors x^^ and x^,. 
can be written as 


r* 

^Isr 


a El 
sr 1 


( 8 - 9 ) 


where a^^ Is a constant which may depend on the spatial distance between the 
r^^ and the s^^ subpixels* Using equations (8-5), (B-6), and (B-9) In 
equation (B-8) yields 



(B-10) 
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The quantity within parentheses In equation (B-IO), In general, depends on the 
number of subpixels and the spatial arrangement of subpixels or on the shape of 
the region of a class In a resolution element. Let £ be a quantity represen- 
tative of the shape of a region of a class In a resolution element: then, 
equation (B-10) can be written as 

» KEj ♦ £KEJ 

« vKiJ (B-11) 

where v « (1 ♦ 6) (B-12) 


If there are R-classes and K.| subcells of each class In a resolution element, 
the spectral vector of the resolution element can be written as 

X - i X.. (B-13) 

1«1 J«1 


The mean vector N of X can be obtained as follows. Consider 

M • E(X) • i E(x. J 
1^ J«1 
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®l”l 


(B-14) 


Assuming the spectral response vectors of subpixels of different classes are 
Independent, the covariance matrix of X can be obtained as follows. 
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Wi 


(B-IS) 


where 84 * r- I 

' '' ( (8-16) 

• (1 + 5^) j 

and 6^ Is a quantity representative of the shape of the reqion of i^^ class in 
a resolution element. A comparison of equations (2-7) and (3-15) shows that 
the effect of correlations between the subpixels of classes is to introduce the 
constants 0 ^. 
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